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ABSTRACT 



Aims. Triton possesses a thin atmosphere, primarily composed of nitrogen, sustained by the sub- 
limation of surface ices 



Methods. The goal is to determine the composition of Triton's atmosphere and to constrain the 



^ ' nature of surface-atmosphere interactions. 

Results. We perform high-resolution spectroscopic observations in the 2.32-2.37 /im range, us- 
ing CRIRES at the VLT. 

Conclusions. From this first spectroscopic detection of Triton's atmosphere in the infrared, we 

, report (i) the first observation of gaseous methane since its discovery in the ultraviolet by Voyager 

OO ■ 

in 1989 and (ii) the first ever detection of gaseous CO in the satellite. The CO atmospheric abun- 
dance is remarkably similar to its surface abundance, and appears to be controlled by a thin, CO- 



(N 
cn 

c7 J ' enriched, surface veneer resulting from seasonal transport and/or atmospheric escape. The CH4 

' partial pressure is several times larger than inferred from Voyager. This confirms that Triton's at- 

mosphere is seasonally variable and is best interpreted by the warming of CH4-rich icy grains as 
Triton passed southern summer solstice in 2000. The presence of CO in Triton's atmosphere also 
affects its temperature, photochemistry and ionospheric composition. An improved upper limit 
on CO in Pluto's atmosphere is also reported. 

Key words. Solar system: general ; Infrared: solar system ; Triton 

1. Introduction 

Like Pluto, Neptune's satellite and probably former Kuiper-Belt object Triton possesses a tenuous, 
predominantly nitrogen atmosphere, in equilibrium with surface ices mostly composed of N2 and 
a variety of other species. The most volatile of these species, CH4 and CO, must be present in trace 
amounts in the atmosphere as well. However, depending on the precise mechanisms of surface- 
atmosphere interactions, the expected atmospheric abundances vary by orders of magnitude, and 
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except for the detection of CH4 in the UV by Voyager in 1989, observations have been severely 
lacking. Progress in IR-detector technology makes the remote study of thin and distant atmospheres 
now possible. Following our observations of methane in Pluto's atmosphere (Lellouch et al. 2009), 
we here report on the first spectroscopic detection of Triton's atmosphere in the infrared. 

2. VLT/CRIRES observations and CH4 and CO measurements 

Spectroscopic observations of Triton were obtained on July 4, 2009, using the CRIRES infrared 
echelle spectrograph (Kaufl et al. 2004) installed on ESO VLT (European Southern Observatory 
Very Large Telescope) UTl (Antu) 8.2 m telescope. We focussed on the regions of the (2-0) band 
of carbon monoxide and of the V3+V4 band of methane, covering the 2318-2330, 2334-2345, 2349- 
2359 and 2363-2373 nm ranges. We used the instrument in adaptive optics mode and with a sUt of 
0.4", providing a mean spectral resolution of 60,000, and acquired spectra during ~ 4 hours. A large 
Doppler shift (-23 km/s) ensured proper separation of the target lines from the telluric absorptions. 

The resulting spectrum shows the detection of many Unes due to methane in Triton's atmo- 
sphere, particularly at 2320-2330 nm (Fig. 1). This is the first observation of gaseous methane 
since its discovery by Voyager (Herbert and Sandel 1991). As for our study of Pluto's CH4, we 
constructed a direct line-by-line atmospheric model of Triton, integrated over angles and including 
solar lines reflected off Triton's surface as well as the telluric transmission (see details in Lellouch 
et al. 2009). The spectrum was first modelled by assuming a single-temperature layer, with Triton's 
atmospheric methane mean temperature (T) and column density (a) as free parameters. We in- 
ferred T=50:^^^ K and a = 0.08±0.03 cm-am (Fig. 2 on-line). The same analysis for Pluto had 
given T=90;|;j3 K and a = 0.75;|;q3q cm-am. This confirms that Pluto's atmosphere is warmer than 
Triton's, as a result of its higher methane abundance. 

The error bars on the inferred mean methane temperature are such that it is not possible to 
constrain the methane vertical distribution. Instead, we used the Voyager-determined thermal struc- 
ture (temperature vs altitude, Krasnopolsky et al., 1993) and methane vertical profile (Herbert and 
Sandel 1991, ingress UV occultation profile). The latter shows a decrease of the CH4 mixing ra- 
tio with altitude with a scale height of ~20 km, due to photolysis. We obtained the same column 
density as above, indicating a partial pressure of methane of 9.8+3.7 nbar, i.e. a surface density of 
(1.9±0.7) X 10'^ cm"-'. This appears to be 4^25 times larger than inferred from Voyager in 1989, 
adopting the CH4 number densities of Herbert and Sandel (1991) and Strobel and Summers (1995) 
(4.7x10'^ cm"-^, within a factor 1.7, averaging ingress and egress). An even larger enhancement 
factor (S'^^) is indicated if the Krasnopolsky and Cruikshank (1995) reanalysis of the Voyager UV 
data, giving CH4 = 3.1±0.8x 10^' cm"^ at the surface, is used. Results are independent on the sur- 
face pressure, as coUisional broadening is negligible. They clearly demonstrate that the CH4 partial 
pressure has increased in the last 20 years. 

The 2335-2365 nm part of the Triton spectrum (see excerpts in Fig. 3) shows the detection of 8 
lines due to the CO(2-0) band (R2-R5, P2, P3, P5 and P8), providing the first detection of CO in its 
atmosphere. An accurate determination of the CO abundance is particularly difficult, as at infinite 
spectral resolution, these features are very narrow, saturated Doppler-shaped lines. Nonetheless, 
assuming a vertically uniform CO distribution, and utilizing the whole set of CO lines (see Fig. 4 
on-line), we determine a CO column of 0.30 cm-am, i.e. a CO partial pressure of 24 nbar, within 
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Fig. 1. Black: Triton spectrum at 2320 - 2330 nm observed by CRIRES/VLT. The spectral resolu- 
tion is 60,000. Red and blue curves show synthetic spectra, including telluric and solar lines. Red: 
methane column density in Triton's atmosphere = 0.08 cm-am. Triton's thermal profile, based on 
Voyager measurements, is taken from Krasnopolsky (1993) and a Voyager-like vertical distribution 
is used for methane (Herbert and Sandel 1991, entrance profile). Blue: no methane. The continuum 
slope is due to the red wing of the + V4 band due to solid methane. Here, as in Fig. 3, the vertical 
unit approximately represents the geometric albedo (but uncorrected for telluric and solar lines). 
Models are shifted vertically by -0.07 for clarity. 

a factor of 3. The column density CO/CH4 ratio is nominally ~3.75 (surface partial pressure ratio 
CO/CH4 ~2.5), with a factor of 4 uncertainty. Deriving the CO/N2 and CH4/N2 mixing ratio is 
complicated by the fact that the surface pressure in 2009 is unknown. Stellar occultation results 
(Olkin et al. 1997, Sicardy et al. 1998, Elliot et al. 1998, 2000a) indicate that the pressure has been 
doubling in ~ 10 years from the 14 yubar value determined by Voyager in 1989 (Gurrola 1995). A 
reasonable assumption for 2009 is 40 //bar, providing CO/N2 ~ 6x10 and CH4/N2 -2.4x10 
at the surface, within factors of 3 and 1 .4 respectively. The CO abundance we determine is many 
times less than previous upper limits (Broadfoot et al. 1989, Young et al. 2001). 

3. Discussion 

Near-infrared observations indicate that CO and CH4 are present on Triton's surface with mixing 
ratios of 0.05 % and 0.1 % relative to N2, and at least for CH4, mostly in solid solution in N2 
(Cruikshank et al. 1993, Quirico et al. 1999, Grundy et al. 2010). In this situation, the expected 
partial pressure of each species is the product of its solid mole fraction and its pure vapor pressure 
(Raoult's law for an ideal mixture). This scenario leads to atmospheric abundances of CO and CH4 
that are about 1 and 3 orders of magnitude lower than observed, respectively (Fig. 5). Although 
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Fig. 3. Black: portions of the Triton spectrum at 2335 - 2340 nm and 2350-2355 nm. Five CO lines 
(R5, R4, R3, P2 and P3) are present in these spectral ranges. Red and blue curves show synthetic 
spectra, including telluric and solar lines, as well as Triton's methane. Red: CO column density 
in Triton's atmosphere = 0.30 cm-am. Blue: No CO. The mismatch in the "continuum" level over 
2351-2354 nm is due to the absorption of the (2-0) band of CO ice (see e.g. Quirico et al. 1999, 
Grundy et al. 2010), not included in models. 



Henry's law may be more applicable than Raoult's in the case of the N2-CH4 system, this does 
not reduce the discrepancy by more than a factor of 2-3. This problem has been studied in the 
case of Pluto's atmospheric methane, present at the 0.5 % level (Young et al. 1997, Lellouch et 
al. 2009), i.e. also considerably enriched over its solid solution equilibrium value. The origin of 
such enhancement is thought to ultimately lie in the seasonal evolution of the Na-dominated solid 
solution. Preferential sublimation of N2 initially creates a thin surface layer enriched in the less 
volatile species. Further evolution of this layer may lead either (i) to the formation of chemically 
pure grains in vertically or geographically segregated deposits ("pure ice" scenario), or (ii) to the 
establishment of a homogeneous "detailed balancing layer" controlling the surface-atmosphere ex- 
changes. In the "pure ice" case, the atmospheric mixing ratios are in simple proportion of the pure 
vapor pressures at the relevant ice temperatures (which may be different for different species) and, 
except for the main species which controls the pressure, of the fractional area covered by each ice. 
Focussing on the case of CH4 on Pluto, Stansberry et al. (1996) demonstrated that pure CH4 lag 
deposits (whose existence is proved by observations, see Doute et al. 1999) assume higher tempera- 
tures than N2 due to their reduced subhmation cooling and preferential formation in regions of high 
insolation. Even if covering only a few percent of Pluto's surface, such patches can explain the ob- 
served atmospheric abundance of methane. Alternatively, the "detailed balanced" model (Trafton, 
1990; Trafton et al. 1998) predicts that surface-atmosphere exchanges in presence of atmospheric 
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Fig. 5. Abundance measurements and equilibrium curves for Triton's volatiles. The black, pink 
and light blue curves show the vapor pressure of pure N2, CO and CH4 ices, calculated from 
Fray and Schmitt (2010). For CO and CH4, the thick curves (red and blue, respectively) show 
the partial pressures based on Raoult's law for an ideal soUd solution with N2, with CO and CH4 
respective abundances of 5 x 10"^ and 1 x 10"^ in the ice (Quirico et al. 1999). For CO and CH4, the 
thin red and blue show the partial pressures expected in the framework of the "detailed balancing 
model" (see text). The black square shows the 1989 Voyager pressure measurement (14 nbar), 
which corresponds to sublimation equilibrium of N2 ice at 37.25 K. The green symbol represents 
the Voyager-measured CH4 partial pressure (Herbert and Sandel 1991); it is plotted at 38+1 K, 
the surface temperature measured by Voyager and the N2 ice temperature inferred by Tryka et al. 
(1994). The blue star and the red square show the CH4 and CO abundances determined in this 
work. They are plotted at a temperature of 39+1 K. A 39 K temperature corresponds to a plausible 
40 fihar pressure for Triton's atmosphere in 2009; 40 K, which corresponds to a N2 vapor pressure 
of 65 yubar, is a reasonable upper limit, given the stellar occultation inferences that Triton's surface 
pressure has doubled in the 10 years following the Voyager epoch. The CO partial pressure we 
measure is consistent with expectations from the detailed balancing model, while CH4 is lower. 



escape and seasonal transport lead to an atmospheric composition reflecting that of the accessible 
ice reservoir from which it is replenished. When no fractionation (e.g. diffusive) occurs during es- 
cape or transport, the atmospheric mixing ratios are identical to those in the volatile reservoir. This 
is accomplished by the thin surface veneer enriched in the less volatile species, throttling off' the 
N2 sublimation, and in permanent equilibrium with the atmosphere according to Raoult's law. 

Our observations provide discriminating keys on these scenarios. The case for CO is most 
straightforward because (i) CO is not subject to diff'usive separation from N2 upon escape (ii) as 
the ratio of its vapor pressure to that of N2 is largely insensitive to temperature in the relevant range 
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Fig. 6. Atmospheric mixing ratios and composition of the ice boundary layer ("film") in the de- 
tailed balancing model. (Top): CO-N2 system. A surface pressure of 40 /ibar is assumed (estimate 
for 2009). The dark blue line is the CO/N2 mixing ratio expected for pure ices. The green curve 
shows the CO mole fraction in the ice surface film, and the black curve is the CO/N2 atmospheric 
mixing ratio derived from the composition of the ice film and applying Raoult's law. The range of 
CO/N2 atmospheric mixing ratios inferred in this work for this pressure ((2-18) x 10 **) is indicated 
by the blue-colored region. It implies a CO/N2 mixing ratio in the surface film of (1.4-12) x 10"^ 
(see text). The surface film is therefore largely dominated by N2, and the total pressure is defined 
by N2 equilibrium at 39.075 K. (Bottom): CH4-N2 system. Calculations are here performed for 
a 14 jubar pressure, appropriate for the Voyager conditions. The Voyager-determined CH4/N2 at- 
mospheric mixing ratio at the surface level ((1.1-3) x 10 **) is indicated by the yellow region. The 
colored lines have the same meaning than in the top panel, with CH4 replacing CO. The intersection 
of the black line with the colored area shows that explaining the observed CH4 mixing ratio and the 
total pressure requires a 38.3-39.6 K surface temperature and would imply a very high CH4 mole 
fraction (50-80 %) in the surface film, well beyond the solubility limit of CH4 in N2. The same 
conclusion is reached if the CH4 amounts measured in 2009 are used. Note that these diagrams 
remain similar at other surface pressures, the only change being the required ice temperature to 
sustain the total pressure. 

(e.g. CO/N2 - 0.112 at 36 K and 0.166 at 41 K), its atmospheric mixing ratio should be roughly 
seasonally constant. The observed CO abundance is over two orders of magnitude lower than the 
pure CO vapor pressure (Fig. 5), and one might envisage that the atmospheric CO results from 
the sublimation of pure CO patches covering ~0.4 % of the surface. However, we do not regard 
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this scenario as likely. Although the pure vs isolated form of CO on Triton's surface has not yet 
been proven from observations, the miscibility of CO and N2 in soUd phase in all proportions and 
the similarity of their vapor pressures argue for a co-condensation of the two species on Triton's 
surface. This is further supported by the strikingly similar longitudinal distribution of the N2 and 
CO ice bands at Triton (Grundy et al. 2010), strongly suggestive of a spatially constant CO/N2 
ice mixing ratio. In addition, even if pure CO patches occurred on Triton's surface, they would 
probably not be able to elevate the CO atmospheric abundance along the mechanism envisaged for 
Pluto methane. This is because CO is not buoyant in N2, restricting the dispersal of CO-rich gas 
in the background atmosphere and therefore the sublimation rates of the CO patches (Stansberry 
et al. 1996). Instead, the detailed balancing model provides a physically-expected interpretation 
to the fact that the atmospheric CO/N2 mixing ratio is consistent with its value in the ice phase 
(Fig. 5). Based on this scenario, the N2-CO composition of the surface boundary layer ("film") can 
be established by simple application of Raoult's law, providing qco(filni) = qco(atm)xp,at(CO) / 
Prat(N2) (Trafton et al. 1998). Adopting again p = 40 //bar, our observed qco(atm) = (2-18)x lO""' 
indicates qco(filni)= (1.4-12)xl0"^. Therefore the surface veneer is still dominated by N2 and the 
presence of CO does not importantly modify the N2 atmospheric pressure, defined by equilibrium 
at 39. 1 K (Fig. 6). Because it may be as thin as a few molecular layers, the surface film may not be 
visible in the near-lR spectra. 

The case for CH4 is more complex. As previously realized (Cruikshank et al. 1993, Yelle et 
al. 1995, Strobel and Summers 1995, Strobel et al. 1996), the CH4 atmospheric mixing ratio at 
the surface measured by Voyager (~1.8xlO"'*) is at least three orders of magnitude larger than 
expected for ideal mixture. However, we note that it is also smaller, by a factor of ~6, than the ice 
CH4/N2 mixing ratio, and as such does not agree with the detailed balancing model in its simplest 
form. Unlike CO, CH4 is subject to atmospheric photolysis and mass separation, and its vapor 
pressure is more temperature-dependent. This probably makes the surface/atmosphere abundance 
relationship for CH4 complex and seasonally variable. In any case, the phase diagram of N2-CH4 
is not obviously consistent with the formation of a CH4-rich solid solution veneer (Stansberry et 
al., 1996). In fact, explaining the range of observed CH4 atmospheric abundance would require 
a CH4 mole fraction in the surface film as high as 50-80 % (Fig. 6), well beyond the solubility 
limit of CH4 in N2 (Prokhvatilov and Yantsevich 1983). The formation of pure CH4 ice grains, 
decoupled from the mixture and not influencing its sublimation (Stansberry et al. 1996, Spencer 
et al. 1997), further evolving into a lag deposit, may be a more plausible outcome. Using a Bond 
albedo of 0.85 (Triton's polar cap) and an emissivity of 0.7-1, a reasonable subsolar temperature 
for these pure methane patches is 45-48 K. Applying the Stansberry et al. (1996) Pluto model then 
indicates that methane patches covering 0.5-1 % of Triton's surface are sufficient to maintain a 
~2xl0"'' atmospheric mixing ratio. Although there is no evidence for such patches in Triton's near- 
IR spectrum, the methane longitudinal distribution of CH4 ice is different from that of N2, and small 
areas of CH4-dominated ice, notably near 300° longitude, are not inconsistent with observations 
(Grundy et al. 2010). In contrast, the existence of widespead pure methane ice is ruled out; therefore 
the emphasized fact that the Voyager-measured methane partial pressure was consistent with vapor 
pressure equilibrium of pure CH4 ice at 38 K is probably coincidental. 

After the Voyager encounter, a variety of seasonal N2 cycle models (see review in Yelle et 
al. 1995) were explored to attempt explaining Triton's visual appearance and the then measured 
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surface pressure. These models, which essentially differed in the assumed ice and substrate albedos 
and thermal inertia, had limited success, leaving unanswered the simple question of where the 
ice is on Triton. Yet, they made distinctive predictions as to the short-term evolution of Triton's 
atmosphere. High thermal inertia models predicted a pressure increase as Triton approached and 
passed Southern summer solstice in 2000 (Triton subsolar latitude moved from 45.5 S in 1989 
to a maximum 50 S in 2000 and 47 S in 2009). This is a consequence of increased insolation 
on, and attendant sublimation of, the Southern polar cap (Spencer and Moore 1992, Forget et al. 
2000). In contrast, "dark frost" models (Hansen and Paige 1992) or low thermal inertia models 
predicted a pressure decrease from ~1980 on, due to the exhaustion of the seasonal southern cap 
and re-condensation of N2 on the invisible winter pole. The discovery of the pressure increase in 
the 1990's, and the persisting signature of N2 and other ices in Triton reflectance spectrum with no 
obvious temporal evolution (Grundy et al. 2010), strongly argue for the fact that the bright deposits 
covering most of Triton southern hemisphere are indeed relatively stable seasonal deposits. Our 
observation that the methane partial pressure has increased by a factor ~4 from 1989 to 2009 is 
qualitatively consistent with the reported pressure increase and the above interpretation. Since the 
CH4 vapor pressure varies more sharply with temperature than N2, we expect that atmospheric 
methane is currently increasing more rapidly than pressure, but multi-volatile seasonal models will 
be needed to fully interpret our results. A direct measurement of Triton's current pressure is also 
highly desirable, and could be obtained through a redetermination of the N2 ice temperature from 
its 2.15 iJ.m band (Tryka et al. 1994). 

The detection of CO has also implications on Triton's atmospheric thermal structure, photo- 
chemistry, and ionosphere. CO is an important cooling agent through radiation in its rotational 
hues (Krasnopolsky et al. 1993, Strobel and Summmers 1995, Elliot et al. 2000b). It eiu^iches 
atmospheric chemistry by introducing additional oxygen species (Krasnopolsky and Cruikshank 
1995). Most importantly, it profoundly modifies ionospheric composition by providing a source of 
C atoms and C"^ ions and by suppressing the N"^ concentration at the benefit of CO"^ and NO"^ (see 
review in Strobel and Sunamers 1995). Although the error bar on the CO abundance is still large, 
all previous considerations on the role of CO have now direct observational support. 

During the same observing night, we also searched for CO in Pluto's atmosphere, covering the 
region of the (3-0) band near 1.57 fim. Only an upper Umit (1 cm-am) was obtained. For a charac- 
teristic surface pressure of 15 jubar (Lellouch et al. 2009), this indicates CO/N2 < 5x10"^. While 

improved over previous results (Bockelee-Morvan et al. 2001, Young et al. 2001), this upper limit 
is still relatively unconstraining when compared to the measured CO ice mole fraction (lxlO~^, 
Doute et al. 1999). Nonetheless, given the similarity of the two bodies, the above considerations 
on the surface control of atmospheric CO at Triton should apply to Pluto as well. In 2015, obser- 
vations with the ALICE and Rex instruments on New Horizons will provide measurements of the 
surface pressure and CH4 and CO abundance in Pluto's atmosphere. We anticipate that CO will 
be measured at a uniform ratio of 0.001 and that the methane mixing ratio will show horizontal 
variability associated with local time and methane patch distribution. 

Acknowledgements. This work is based on observations performed at the European Southern Observatory (ESO), proposal 
383.C-0609. We thank F. Forget, B. Schmitt, and L. Trafton for illuminating discussions. 
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Fig. 2. Zoom on several CH4 lines showing sensitivity of the spectrum to Triton's methane. Some 
portions of the spectrum, especially over 2320-2326 nm, are relatively independent on tempera- 
ture, while high-energy lines at 2353-2359 nm show increased temperature sensitivity.The top two 
panels show sensitivity to the methane abundance. Blue, red, and green synthetic spectra have 0.03, 
0.08, and 0.20 cm-am of methane. Triton's thermal profile is taken from Krasnopolsky (1993) and 
a Voyager- like vertical distribution is used for methane (Herbert and Sandel 1991, entrance pro- 
file). Based on these models, the best fit methane column density is determined to be 0.08+0.03 
cm-am (i.e. +40 %). The bottom two panels show sensitivity to methane temperature. The previous 
best-fit model using Voyager thermal profile is shown in pink. Other models assume an isothermal 
atmosphere with temperature of 30 K (dark blue), 50 K (red) and 80 K (green). These fits indicate 
a mean methane temperature of 50^j5 K. 

Figure available on-line. 



2358.8 



E. Lellouch et al.: Detection of CO in Triton's atmosphere, Online Material p 2 



1.1 
1 

0.9 
0.8 
0.7 
0.6 



1 


1 


R5 


No CO 




0.03 cm-am 




0.3 cm-am 


1 


3 cm-am 

1 





2335.1 2335.3 2335.5 2335.7 



2337 2337.2 2337.4 2337.6 



2338.8 2339 2339.2 2339.4 




2340.6 2340.8 2341 2341.2 




2357.9 2358.1 2358.3 2358.5 




2351.2 2351.4 2351.6 2351.8 
1.05 



1 1 



0.85 
0.8 
0.75 



P8 



2365 2365.2 2365.4 2365.6 




2353.4 2353.6 2353.8 2354 



Fig. 4. The eight detected CO Hnes from Triton's atmosphere. X-axis units are nm and Y-axis units 
are arbitrary. Lines are compared with models including 0, 0.03, 0.30 and 3 cm-am of CO. The 
slow change of absorption depth with change of column density is caused by heavy saturation 
of particularly narrow Doppler-shaped lines at T ~50 K. Based on these models, the best fit CO 
column density is determined to be 0.30 cm-am with a factor of 3 uncertainty. 

Figure available on-line. 



